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ABSTRACT 

A  computer  program  using  a  curve-fit  technique  was  developed  to  reduce 
performance  data  obtained  from  static  tests  of  aircraft  propellers.  The  entire 
program  is  written  in  Fortran  IV  language  for  use  on  the  IBM  7094  computer 
located  at  Wright-  Patters  on  AFB,  Ohio. 

The  program  accepts  static  whirl  rig  test  data  (ie,  raw  RPM,  horsepower, 
and  thrust  data)  obtained  at  a  fixed  blade  angle  and  reduces  it  into  pertinent 
propeller  relationships.  The  program  first  reduces  the  test  data  into  various 
coefficients  and  computes  the  propeller  tip  Mach  number.  A  curve  fit  technique 
hen  fits  running  curves  through  the  test  thrust  and  horsepower  data  points  at  the 
test  tip  Mach  numbers.  Intermediate  horsepower  and  thrust  values  are  determined 
from  the  fitted  curves  at  selected  Mach  number  increments,  and  all  coefficients 
a  re  recomputed.  This  results  in  a  presentation  of  the  reduced  data  in  two  forms, 
ooeffclents  computed  from  the  actual  test  data  and  coefficients  obtained  at  specific 
constant  Mach  number  increments  from  the  fitted  curves.  The  data  is  presented 
in  tabular  printout  form.  This  is  a  general  program  and  is  written  so  that  the 
order  of  the  curve  fit,  the  Mach  number  increment,  the  number  of  test  data  points 
and  the  Mach  number  range  can  be  varied. 

The  entire  program  deck  and  all  nonstandard  subroutines  are  included. 
Detail  instructions  are  provided  which  should  allow  the  program  to  be  used  by 
technicians  or  students  who  are  familiar  with  Fortran  IV  language. 
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SECTION  I 
INTRODUCTION 

This  program  was  developed  to  reduce  data  accumulated  during  an  extensive 
series  of  static  tests  of  propellers  for  the  XC-142A  Y/STOL  aircraft  conducted 
at  Wright- Patterson  Air  Force  Base,  Chlo.  These  tests  are  described  in  detail 
in  a  report  being  prepared  entitled  “Propeller  Static  Performance  Tests  for 
V/STOL  Aircraft.”  The  program  is  complete  in  that  it  contains  all  non-standard 
subroutines. 

The  program  reduces  RPM  and  corrected  horsepower  and  thrust  data  at 
any  given  blade  angle  and  computes  Power  Coefficient  (Cp),  Thrust  Coefficient 
(Ct),  C^Cp,  Figure  of  Merit  (F.M.)  ani  thrust/horsepower  (TH/HP).  Propeller 
tip  Mach  number  is  calculated  using  propeller  diameter,  RPM  and  ambient 
temperature.  A  curve  fit  technique  is  then  employed  to  fit  a  running  selected 
degree  polynomial  using  least  squares  to  NPHP  (horsepower)  or  NPTH  (thrust) 
consecutive  test  data  points  (see  NAMELIST)  at  the  test  tip  Mach  numbers.  In 
this  case,  a  2nd  degree  polynomial  was  fit  to  6  consecutive  data  points.  Inter¬ 
mediate  horsepower  and  thrust  values  are  determined  from  the  fitted  curves  at 
the  selected  Mach  number  Increments.  The  routine  proceeds  until  smoothed 
curves  for  horsepower  (HPI)  and  thrust  (THI),  and  derived  coefficient  values  plus 
corresponding  RPM  and  tip  speed  (TIPS),  have  been  computed  at  all  of  the  selected 
increments  of  Mach  number.  The  procedure  is  essentially  one  of  creating  equal 
tabular  entries  through  a  smoothing  technique  rather  than  through  an  lnter- 
polative  procedure. 

The  Namelist  variables  establish  the  boundaries  of  the  program,  such  as 
the  Mach  number  increments  and  range,  number  of  data  points  fitted  by  the 
polynomial,  and  the  degree  of  the  polynomial.  All  variables  can  be  changed  to 
adapt  the  program  to  other  uses. 
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section  n 

COMPUTER  PROGRAM 

1.  PROGRAM  OPERATION 

The  computer  program  is  designed  to  provide  parameters  for  aircraft 
propellers  at  specified  tip  Mach  number  increments  by  fitting  curves  to  data 
for  horespower  and  thrust  versus  Mach  number  at  a  set  blade  angle  and  re¬ 
computing  values  for  the  parameters  at  the  predetermined  Mach  number  in¬ 
crements.  A  flow  chart  showing  the  order  of  the  program  and  the  data  input 
is  given  in  Figure  1.  Data  is  obtained  from  a  test  rig  data  sheet,  such  as  is 
shown  in  Figure  2,  and  coded  for  the  computer  as  shown  in  Figure  3. 

For  a  given  blade  angle,  the  data  must  be  presented  to  the  compute*,  at 
increasing  values  of  RPM.  After  receiving  the  data  input,  the  computer  cal¬ 
culates  the  various  parameters  from  the  test  data,  including  Mach  number,  C^, 
Cp,  C/C  .  F.M.,  and  TH/HP,  according  to  the  appropriate  equations  (given  in 
the  list  of  program  variables).  Values  are  printed  out  on  the  output  sheet. 
Figure  4,  as  Raw  Data  Points. 

Subroutine  XIPLSQ  is  then  called  to  fit  a  6-point,  2nd  order,  least  squares 
curve  through  the  data  points  for  HP  and  TH  at  test  tip  Mach  numbers.  Then  the 
computer  determines  values  for  HP  and  TH  from  the  fitted  curves  that  cor¬ 
respond  to  the  selected  Mach  numbers,  and  uses  these  values  to  recompute 
values  for  all  parameters  at  the  selected  Mach  numbers.  These  values  are  then 
printed  out  on  the  output  sheet  as  shown  on  Figure  4  under  Fitted  Curve  Data 
for  Constant  Mach  Number  Increments. 

The  computer  then  reads  in  the  data  for  a  new  blade  angle  (2nd  Beta, 
3rd  Beta,  etc.)  and  repeats  the  computations. 

2.  SUBROUTINE  XIPLSQ 

The  least  squares  curve  is  fitted  by  means  of  Subroutine  XIPLSQ.  HPI  and 
THI  values  are  computed  as  functions  of  the  independent  variable  BMACH  by 
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ORDER  OF  INPUT  SEQUENCE 

Cord  Column  I 


PROGRAM  DECK 


1st 

PROP' 


1st  BETA 


2nd  BETA 


LAST  BETA 


LAST_ 

prop" 


1st  BETA  - 


LAST  BETA 


$  JOB 


END 


XIPLSQ  Subroutine 
(DKEND  XYPLSQ 
$DATA 
$NAM  $ 


PROPELLER  IDENTIFICATION  CARD 

PROPELLER  AND  ATMOSPHERIC  DATA 
10  20 


1st  RPM 


HP 


Last  RPM  HP 

BLANK 


IDENTIFICATION  CARD 

Last  RPM  HP 

BLANK 


IDENTIFICATION  CARD 


Last  RPM 
-I . 


HP 


IDENTIFICATION  CARD 


Lost  RPM  HP 

BLANK 


IDENTIFICATION  CARD 


Last  RPM 
- 1 , 


HP 


END  OF  FILE  CARD 


30 

TH 


TH 


TH 


TH 


TH 


TH 


Figure  1.  Sample  Input  Data 
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Figure  2.  Typical  Rig  Data  Sheet 


FORTRAN  COOING  FORM 


ASD-TR-68-l'J 


Figure  3.  Sample  Coding  Form 

(Taken  from  Typical  Rig  Data  Sheet) 


>FH  li  Al-4  A  7J  V-IJ6  'tlH  .0  230  Vi-ULS  POWM  ie>000  PT-LP 
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Figure  4.  Typical  Computer  Printout 
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means  of  a  polynomial  smoothing  process;  this  process  resembles  nonlinear 
interpolation  in  that  it  considers  several  points  of  the  HP  and  TH  versus 
A  MACH  tabular  data  function  on  either  side  of  the  desired  BM4CH  value  in 
determining  a  least  squares  polynomial,  and  this  polynomial  is  then  used  to 
determine  the  smoothed  HPI  and  THI  values. 

Smoothing  in  the  vicinity  of  the  ends  of  the  test  data  array  will  not  produce 
as  satisfactory  results  as  in  the  center  because  there  are  not  sufficient  data 
points  available  on  both  sides  of  the  desired  BMACH  argument.  The  polynomial 
is  always  fit  to  the  predetermined  number  of  test  data  points,  however,  so  addi¬ 
tional  data  points  from  the  other  side  are  used  to  augment  the  least  squares 
input  data.  If  the  desired  BMACH  value  is  out  of  range  of  the  test  data,  however, 
extrapolative  smoothing  would  be  required,  which  would  give  results  of  even 
less  reliability. 

Note  that  if  the  fitted  polynomial  is  of  the  Nth  degree,  N  +  1  data  points  are 
required  to  develop  the  smoothing  polynomial. 

3.  DATA  INPUT 

The  first  data  input  card  (propeller  identification  card)  shows  which  propeller 
data  is  being  computed.  The  second  card  gives  the  basic  propeller  data  and 
atmospheric  conditions  at  the  time  that  the  test  was  run,  including  the  number  of 
blades,  propeller  diameter,  activity  factor,  blade  angle  setting,  temperature  (°C), 
and  density  ratio  (air  factor).  The  first  RPM  value  with  its  corresponding  HP 
and  TH  values  is  punched  on  the  third  card,  and  values  for  increasing  RPMs  are 
punched  on  succeeding  cards.  A  blank  card  is  placed  after  the  last  RPM-HP-TH 
card  if  another  blade  angle  run  is  to  be  computed  for  that  propeller;  if  no  more 
blade  angle  runs  are  to  be  made,  a  “-1.”  card  (with  -1.  punched  within  the  first 
10  spaces)  follows. 
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SECTION  m 

PROGRAM  PARAMETERS 

1.  NAMELIST 

The  variables  in  the  Namelist  control  various  aspects  of  the  program.  They 
are  set  at  specified  values  by  the  DATA  statement  near  the  beginning  of  the 
program  deck. 

There  is  one  Namelist,  NAM,  for  this  program,  consisting  of  the  following 
variables: 

VL  =  VL  +  VS,  -  minimum  selected  Mach  number  for  curve- 
fitted  data  (VL  set  at  0.500) 

VH  -  maximum  selected  Mach  number  for  curve  fitted  data 
(VH  set  at  1.000) 

VS  -  increment  increase  in  selected  Mach  numbers 
(VS  set  at  0.025) 

NPTH  -  maximum  number  of  raw  thrust  points  considered  at 

one  time  for  a  curve  fit  (NPHP  set  at  6) 

NPHP  -  maximum  number  of  raw  horsepower  points  considered 

at  one  time  for  a  curve  fit  (NPHP  set  at  6) 

NDTH  -  order  of  the  polynomial  used  for  fitting  the  thrust  curve. 

(NDTH  is  set  at  2;  that  is,  a  running  polynomial  of  the 

2 

form  CQ  +  Cj  x  +  C2  x  is  used  in  smoothing  the  thrust 
data  unless  modified  by  parameter  control.) 

NDHP  -  see  information  for  NDTH,  above  (set  at  2). 

Namelist  values  can  be  modified  by  the  simple  procedure  of  replacing  one 
card.  For  example,  if  the  value  for  VH  (maximum  Mach  number)  is  to  be  changed 
from  1.000  to  0.900,  the  following  card  would  be  removed 

Card  Column  2 

$  NAM  $  (the  first  card  after  the  $  DATA  control  card) 

and  the  following  card  inserted: 

Card  Column  2 

$  NAM  VH  =  .9$ 
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2.  EQUATIONS 

The  following  equations  are  used  in  computing  the  parameters. 


TR  =  1.  8(TC  +  273. 16) 


AMACH  '  (49™!  ^ 


CT  =  CT(TH)  =  0.  1518  x  10 


TH 


(RPM)2  (DIA)4 


CTI  =  CT(THI) 


(see  THI) 


CP  =  CP(HP)  =  0.5  x  10 


CPI  =  CP(HPI) 


11 


HP 


(RPM)3  (DIA)5 


(see  HPI) 


ICTl 


1.5 


FM  rr  FM(CT,  CP)  =  0.  798  —  ;  if  CT  <  0,  FM  is  set  at 

value. 


FMI  FM(CTI,  CPI) 


3.  OTHE  R  PROGRAM  VARIABLES 


a.  Nonsubscripted  Variables 


NBL  - 
AF  - 
DIA  - 
RPM  - 
CT  - 
CP  - 
FM  - 
TIPS  - 


number  of  blades  per  propeller  (reference  data) 

activity  factor  (reference  data) 

diameter  of  propeller,  feet 

propeller  speed  (revolutions  per  minute) 

coefficient  of  thrust 

coefficient  of  power 

Figure  of  Merit 

actual  tip  speed,  ft/ sec 


minus 
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b.  Subscripted  Variables 


HEAD 

B 

TC 

S 

TR 

HP 

TH 

AMACH 

BMACH 

THI 

HPI 

FMI 

CTI 

CPI 


propeller  identification 

BETA,  blade  angle,  degrees  (reference  data) 

TEMPC,  ambient  temperature  (°C) 

"sigma"  =  (air  factor)  density  ratio  (a) 

TEMPR,  ambient  temperature  (°R) 

corrected  horsepower;  (Corrected  HP  =  — ) 

Test  TH 

corrected  thrust,  pounds  (Corrected  TH  =  — — — ) 

actual  tip  Mach  number 

selected  Mach  number  for  curve  fit 

thrust  values  determined  irom  the  smoothed  thrust  curve. 

horsepower  values  determined  from  the  smoothed  horse¬ 
power  curve. 

interpolated  figure  of  merit 
interpolated  coefficient  of  thrust 
interpolated  coefficient  of  power 


c.  Counters 


NP 

NB 

N 

NV=J-1 


number  of  propellers 

number  of  blade  angles 

number  of  cards  per  blade  angle  run 

number  of  selected  Mach  numbers  per  run 


d.  Miscellaneous  Parameters 


AF 


activity  factor  = 


100,000 

f  b  1 

f 

3 

I  d  1 

fJL.  \ 

16  J 

D 1  A  * 

v  R  1 

{  R  / 

0.2 
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b  -  blade  width,  feet 

R  -  total  blade  radius,  feet 

r  -  radius  along  blade,  feet 

<y  -  density  ratio;  a  =  — —  (reference  data) 

Po 

p  -  local  density,  lb  sec  /ft4 

p  -  sea  level  standard  density,  lb  sec  /it* 
o 
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SECTION  IV 

COMPARISON  OF  CURVES 

A  comparison  of  the  computer  curve-fit  data  with  hand-faired  data  for  thrust 
and  horsepower  is  shown  in  Figures  5  and  6.  These  figures  show  that  the  test 
data  and  the  computer  curve-fit  provide  very  good  agreement.  The  curve-fit 
results  in  a  substantial  smoothing  of  the  data. 


12 
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Figure  5.  Comparison  of  Thrust  Curve- Fit  Data  to  Raw  Data 
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MACH  NUMBER 


Figure  6.  Comparison  of  Horsepower  Curve-Fit  Data  to  Raw  Data 
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APPENDIX 

COMPUTER  PROGRAM 

This  appendix  presents  the  basic  computer  program 
and  the  necessary  subroutines,  followed  by  an  expla¬ 
nation  of  the  subroutines. 
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EXPLANATION  OF 
SUBROUTINE  MTXEQ 

MTXEQ,  Matrix  Equation  Solver,  solves  the  system  of  linear  equations 
expressed  by  the  matrix  equation  AX  =  B,  using  single  precision  floating 
arithmetic. 

Control: 

To  solve  the  matrix  equation, 

AX  =  B 

the  following  calling  sequence  is  used: 

DIMENSION  A(N,N),  X(N,K),  B(N,K) 

CALL  MTXEQ(A.X,B,N,K) 

where 

A  -  is  the  N  by  N  coefficient  matrix.  The  A  matrix  is  not 
destroyed. 

X  -  will  be  the  N  by  K  solution  matrix. 

B  -  is  the  N  by  K  right-hand-side  matrix.  The  B  matrix  is 
not  destroyed. 

N  -  is  the  number  of  equations,  i.e.,  rows  in  A,  X  and  B. 

(1  <  N  <  24) 

K  -  is  the  number  of  solution,  i.e.,  columns  in  X  and  B. 

(K  >  0  and  (N+K)  <  26) 

Other  Programming  Information: 

1.  Other  subprograms  used  -  None. 
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2.  Error  conditions  - 

a.  The  following  ranges  of  N  and  K  are  allowed: 

(1  <  N  <  24) 

K  >  0  and  (N+K)  <  26 

When  violations  of  the  above  ranges  are  detected  by  MTXEQ, 
a  message  is  written  and  FXEM  is  called.  The  upper  limits 
on  N  and  (N+K)  may  be  altered  by  changing  a  DIMENSION 
and  a  DATA  statement. 

b.  If  the  A  matrix  has  DET(A)  =  0,  a  message  is  written 
and  FXEM  is  called. 

c.  Common  usage  - 

MTXEQ  uses  664  cells  of  blank  common. 

Method: 


The  matrix  C  =  [A,Bj  is  formed  in  blank  common.  Subsequent  operations 
are  performed  on  the  C  matrix,  thus  preserving  the  A  and  B  matrices. 
Elementary  row  operations  are  performed  to  eliminate  elements  below  the 
main  diagonal  of  the  augmented  matrix  C.  To  reduce  the  propagation  of  round¬ 
off  error,  the  i-th  pivot  element  is  chosen  to  be  an  element  having  maximum 
magnitude  in  the  i-th  sub-column.  Finally,  a  back  solution  gives  the  desired 
matrix,  X.  The  number  of  floating  point  arithmetic  operations  required  using 
this  method  is  of  the  order  of  2/3  (N2  +  3KN2). 
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EXPLANATION  OF 
SUBROUTINE  PLSQ 


PLSQ,  Polynomial  Least  Square  Curve  Fit,  will  fit  a  given  set  of  data  to  a 
polynomial  of  degree  K  of  the  form: 

Y  =  C1Xk+C2Xk“1+...+CkX+Ck+1 

An  error  analysis  of  the  fit  is  optionally  given. 


Control: 


DIMENSION  X(N),  Y(N),  C(L) 

where  L  is  K+l 

CALL  PLSQ(X,Y,N, K,C,LIST,EMAX,ERMS,EMEQ) 

where 

X  is  the  array  of  N  independent  variables. 

Y  is  the  array  of  dependent  variables. 

N  is  the  number  of  independent  (dependent)  variables. 

K  is  the  degree  of  the  least  squares  polynomial. 

C  is  the  array  of  the  coefficients,  high  order  to  low 
order,  of  the  least  squares  polynomial. 

LIST  =  0  suppresses  the  error  analysis  output. 

=  1  gives  the  error  analysis  output. 

EMAX  is  the  maximum  absolute  error  obtained  by  using  the 
computed  C  coefficients  to  approximate  the  dependent 
variable. 

ERMS  is  the  root  mean  square  error  obtained  by  using  the 
computed  C  coefficients  to  approximate  the  dependent 
variable. 

EMEQ  is  the  maximum  deviation  from  unity  in  the  linear  system 
check  solution. 
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(Xher  Programming  Information: 


1.  Other  subroutines  used: 

PLSQ  calls  MTXEQ  to  solve  the  normal  equations  and 
FXEM  for  return  to  the  system  on  error  conditions. 

2.  Error  conditions: 

PLSQ  checks  for  the  following  errors: 

a.  K  >  23 

b.  N  <  K 

c.  K  <  O 

Upon  detecting  an  error  a  pertinent  message  is  given 
and  FXEM  is  called  for  error  tracing  and  return  to 
the  system. 

3.  Common  usage: 

PLSQ  uses  1309  cells  of  blank  common. 

Method: 


PLSQ  fits  a  set  of  observed  data,  [Xj»yjJ»  with  a  polynomial  of  the  form: 
Y  =  C1Xk+C2Xk"1+...+CkX+Ck+1 

by  solving  the  normal  equations  for  the  C.’s.  A  transformation  of  the  range  of 
X  to  (-1,1)  gives  a  partial  decoupling  of  the  normal  equations  which  improves 
the  accuracy  of  the  solution.  Double  precision  arithmetic  is  performed  at 
critical  points  in  the  computation. 
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112.  SPONSORING  MILITARY  ACTIVITY 


Deputy  for  Engineering 

Wright- Patterson  Air  Force  Base,  Ohio 


A  computer  program  using  a  curve-fit  technique  was  developed  to  reduce 
performance  data  obtained  from  static  tests  of  aircraft  propellers.  The  entire 
program  is  written  in  Fortran  IV  language  for  use  on  the  IBM  7094  computer 
located  at  Wright- Patterson  AFB,  Ohio. 

The  program  accepts  static  whirl  rig  test  data  (ie,  raw  RPM,  horsepower, 
and  thrust  data)  obtained  at  a  fixed  blade  angle  and  reduces  it  into  pertinent 
propeller  relationships.  The  program  first  reduces  the  test  data  into  various 
coefficients  and  computes  the  propeller  tip  Mach  number.  A  curve  fit  technique 
then  fits  running  curves  through  the  test  thrust  and  horsepower  data  points  at  the 
test  tip  Mach  numbers.  Intermediate  horsepower  and  thrust  values  are  determined 
from  the  fitted  curves  at  selected  Mach  number  increments,  and  all  coefficients 
are  recomputed.  This  results  in  a  presentation  of  the  reduced  data  in  two  forms, 
coeffcients  computed  from  the  actual  test  data  and  coefficients  obtained  at  specific 
constant  Mach  number  increments  from  the  fitted  curves.  The  data  is  presented 
in  tabular  printout  form.  This  is  a  general  program  and  is  written  so  that  the 
order  of  the  curve  fit,  the  Mach  number  increment,  the  number  of  test  data  points 
and  the  Mach  number  range  can  be  varied. 

The  entire  program  deck  and  all  nonstandard  subroutines  are  included. 
Detail  instructions  are  provided  which  should  allow  the  program  to  be  used  by 
technicians  or  students  who  are  familiar  with  Fortran  IV  language. 
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